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ABSTRACT 

Context. Studying Luminous InfraRed Galaxies (LIRGs) is particularly important in the build-up of the stellar mass from z=l to 
z=0, and to determine physical properties of these objects at the redshift 0.7. LIRGs are now identified to play a major role in galaxy 
evolution from z=l to 0. The global star formation rate (SFR) at z ~ 0.7 is mainly produced by LIRGs. 

Aims. We perform a multiwavelengths study of a LIRGs sample in the Extended Chandra Deep Field South at z=0.7, selected at 
24 jim by MIPS onboard S pitzer S pace Telescope and detected in 17 filters. Data go from the near-ultraviolet to the mid-infrared. 
This multiwavelengths dataset allows us to bring strong constraints on the spectral energy distributions (SEDs) of galaxies, and thus 
to efficiently derive physical parameters as the SFR, the total infrared luminosity, attenuation parameters and star formation history. 
We distinguish a sub- sample of galaxies detected at 70 iim which we compare to the rest of the sample to investigate the relative 
importance of this wavelength on the determination of the physical parameters. An important part of this work is the elaboration of a 
mock catalogue which allows us to have a reliability criteria for the derived parameters. 

Methods. We study LIRGs by means of a SED-fitting code CIGALE. At first, this code creates synthetic spectra from the Maraston 
(2005) stellar population models. The stellar population spectra are being attenuated by using a synthetic Calzetti-based attenuation 
law before the addition of the dust emission as given by the infrared SED library of Dale&Helou (2002). The originality of CIGALE 
is that it allows us to perform consistent fits of the dust-affected ultaviolet-to-infrared wavelength range. This technique appears to be 
a very powerful tool in the case where we can have access to a dataset well-sampled over a large range of wavelengths. 
Results. We are able to derive a star formation history and to estimate the fraction of infrared luminosity reprocessed by an active 
galactic nucleus. We study the dust temperatures of our galaxies detected at 70 fim and find them colder than predicted by models. 
We also study the relation between the SFR and the stellar mass and do not find a tight correlation between both, but a flat distribution 
and a large scatter which is interpreted in terms of variations of star formation history. 
Conclusions. 
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1. Introduction 

Luminous InfraRed Galaxies (LIRGs) are commonly defined 
as galaxies whose infrared (IR, 8-1000 jim) emission is higher 
than 10 n L o and lower than 1O 12 L . In a pionneering work 



|Chary & Elbaz] ( [2001) (hereafter CE01), identified them as the 
main contributor to the cosmic IR background detected by IS O 
(Infra red S pace Observatory) ; the deep S pitzer/MlPS (Rieke 
|et al. | ( [2004] )) surveys have confirmed that they play a major role 
in galaxy evolution at intermediate redshift (z=0.5-l). Indeed the 
IR luminosity function is found to evolve with z and to shift to- 
wards higher luminosities when z increases (Le Floc'h et al. 
2005), |Pe7e^ Gonzale z J eTai71 ([2005)) as a consequence, LIRGs 



are rare objects in the nearby universe but become the major con- 
tributors to the star form a tion density measured in IR at z=0.5-l 
Le Floc'h eTaT~|p005] ), |Caputi et ai~|p007] ), [Magnelli et al" 
2009 )). The evolution of LIRGs as a function of z does not only 



concern their contribution to the luminosity density but also their 
physical characteristics. Most of the LIRGs populating the lo- 
cal u niverse are located in systems interacting violently (Ishida 
( [2004] )). Conversely, at z ~1 only 30% of LIRGs exhibit fea- 



tures linked to violent merging (Bell et al. ( 2007 ), Zheng et al. 
(2007)) : most of them look like bright spirals who experience a 
secular evolution without violent events. This morphological dif- 
ference between local and distant LIRGs is corroborated by the 
analysis of their star formation rate (SFR). Whereas local LIRGs 
are experiencing a strong starburst, distant LIRGs do not seem 
to strongly depart from the mean SFR - stella r mass (hereafter 
Af*) relation found at z=l ( |Elbaz et al.] ( |2007] )). 

Quite naturally, most of the studies devoted to LIRGs are 
based on IR data only. Very recently several teams aimed at 
studying local LIR Gs by com bining a larg e set of w avelengths 
and observations Burgarel la et al. ( 2005| ), Kaviraj et al. (2008), 
|Da Cunha~| ( |2008| )). At higher z such studies are scarcer. Zheng 
et al. (2007) have combined ultraviolet (UV), optical and IR ob- 
servations and built mean spectral energy distribution (SEDs) of 
IR selected galaxies at z=0.7. They are found similar to those 
of nearby normal galaxies with lower IR luminosities (hereafter 
L dusu L(8 - 1000/mi)). 

The availability of large datasets of photometric data have 
lead to a renewal of stellar populations synthesis codes assuming 
more or less complex star formation histories (SFHs). The com- 
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bination of galaxy simulations with realistic SED modelling im- 
proves our understanding of uncertainties of these fitting meth- 
ods to estimate SFRs, M+ or population ages dWuyts et al. 
( [2QQ^ , |Conroy et al.| ( |2Q09a] , |Lee et al.lpOD^ , |Dye | ( |20Q8| ))7 
Anyway, most of these studies only deal with the stellar UV - 
Near-InfraRed (NIR) SEDs of galaxies. Since LIRGs emit most 
of their energy in IR we must account for the entire UV-to-IR 
SEDs to characterise them. Up to now there has been few at- 
tempts to combine stellar and dust emissions in the modelling 
of galaxy SEDs. Stars and dust interact in a very complex way 
in galaxies and a detailed description of this interplay is not 
always consistent with a SED-fitting code aim ed at studying 
large samples of galaxies (Dopi ta et al. (|2005|)). The GRASIL 



January as part of the MIPS Guaranteed Time Observing pro- 
grammeme. 



(GRAphite and SILicate) code (Su va et al. | f!998), P anuzzo et| 



|al. | ( [2005] ) ) produces UV-to-IR SEDs ba sed on a set of physical 
processes; Iglesias-Paramo etai7] ( |2QQ7| ) used GRASIL to gener- 
ate a synthetic library of galaxy SEDs and fit a multiwavelengths 
dataset of UV selected galaxies from z=0.2 to 0.7. Da Cunha et 
al. (2008) developed a physically motivated model combining 
three IR components to the stellar populations synthesis code 
of|Bruzual & Chariot |([2003]). Burgarella et al. (2005) and [Noll 



et al. '( 2009b| ) followed a similar approach and proposed a code 



which calculates the effect of dust on galaxy SEDs in a consistent 
way; the code is versatile enough to be used at different redshifts 
and with different SFHs. All these codes which combine stellar 
and dust populations are very efficient to constrain dust attenu- 
ation on the basis of an energetic budget : the IR dust emission 
originates from the stellar light absorbed in the UV-NIR part of 
the spectrum. In this paper we apply a multiwavelengths anal- 
ysis from the far-ultraviolet (FUV) to the IR, based on SED- 
fitting, on a sample of z=0.7 LIRGs selected at 24 yum. Our anal- 
ysis is based on the code CIGALE (Code Investigating GALaxy 
Emission^} ( [Noll et a~ ( 2009b )). It is the first application of 
this code to intermediate redshift galaxies mostly emitting in IR. 
Our aim is to study this galaxy sample, representative of LIRGs 
at intermediate redshift, in a very homogeneous and systematic 
way to determine the main characteristics of their stellar popula- 
tions and dust emission. We will re-visit dependance of the dust 
temperature (deduced from the IR SED) on L^ust and the SFH of 
these galaxies. 

In section 2 we describe the sample and the multiwave- 
lengths dataset. In section 3 we summarize the main character- 
istics of the CIGALE code and its application to our sample, 
in particular we build a mock catalogue to test the code perfor- 
mances. The results of the SED-fitting to the LIRGs sample are 
presented in section 4. Section 5 is devoted to a discussion of 
these results in the framework of the temperature-luminosity re- 
lation and the specific SFRs of LIRGs at z=0.7 which will be 
compared to model predictions. 

We adopt an H () =70 km s" 1 Mpc" 1 , Q M =0.3, Q A = 0.7 cos- 
mology throughout this paper. 



2. The sample and data 



The C handra Deep Field South (CDFS, Giaccon i et al. 
(2002)) is one of the most intensely studied regions of the 
sky, with observations stretching from the X-ray to the ra- 
dio, including UV, optical, IR, and submillimeter imaging, 
from space-based as well as the largest terrestrial observato- 
ries. The region of the CDFS (<*=3h32m00s, £=-27° 35' 00", 
J2000.0) was observed with the MIPS instrument on board 
the S pitzer S pace Telescope (Wern eFet al. | ( |2004| )) in 2004 



|Le Floc'h et al. | ( |2005| ) used the 2955 sources detected at 24 
fim by MI PS with flux at 24 /im > 83 /iJy at the 80% complete- 
ness limit (Pap ovich et al. | ( |2004| )) to build the IR luminosity 
function from z=l to 0. Photometric redshifts from COMBO- 17 
(Classifying Objects by Medium-Band Observations in 17 fil- 
ters, [WoireFal7] ([2004]) ) were used for sources at z < 1.2 (see 
Le Floc'h et al. (2005) for more details). [Buat et al.| ( |2007a| ) se- 
lected a subsample of 190 LIRGs at 0.6<z<0.8 and measured 
the near ultraviolet (NUV) emission of these galaxies on deep 
Galaxy Evolution Explorer (GALEX) images. 80% of the 190 
LIRGs were detected at NUV. The four IRAC bands (3.6/mi, 
4.5yum, 5.8/mi and 8.0/mi) are also available for these sources 
(N. Bavouzet, private communication). The optical-to-NIR pho- 
tometry was taken from a K- selected catalogue of the ECDFS( 
Extended CDFS, Of=3h32m00s, £=-27° 48' 00") a s part of the 
multiwavelengths Survey by Yale-Chile (MUSYC, |Taylor et al. 
fl([2OO9)0). This catalogue consists of photometry derived from 
U, U38, B, V, R, I, z', J, H and K, imaging for approximately 
80% of the field. The 5 <x flux limit for point-sources is 
=22.0. 



We cross-correlate sources selected by |Buat et aLj p007a| ) 
and detected at 24 /mi with MUSYC data using a tolerance ra- 
dius of 2". When a double match is found, we select the closest 
object. We ignore the case where three or more optical sources 
can be associated with a given 24 /mi detection and thus elimi- 
nate 9 objects from the 190 sources. At the end we are left with 
181 galaxies whose 150 have a NUV detection (Table [T]) with a 
mean redshift of 0.70 +/-0.05. For the 31 objects not detected in 
NUV a mean redshift of 0.72 +/- 0.05 is found. 

The ECDFS was observed with the MIPS instrument as part 
of the FIDEL (Far-Infrared Deep Extragalactic Legacy Survey, 
PI : Dickinson) legacy programmeme . For the 70 jum data, we 
used the catalogue from Magne lTTet al. | ( [2009| ) . Detection limit 
for 70 fim sources in this field is 3.5 mJy (5 <x). We refer to 
this paper for more details about the source extraction, photom- 
etry and completeness of the catalogue. We cross-co rrelated our 
181 objects with FIDEL data thanks to their IRAC ( Fazio et al. 
[] ( |2004| )) positions, with a a tolerance radius of 2". Finally 62 
galaxies were detected at 70 yum (Table [T}, with a mean redshift 
of 0.68 +/- 0.04. 

The fraction of MIPS sources showing evidence for the pres- 
ence of an Active Galactic Nucleus (AGN) in their optical coun- 
terparts is less than -15% according to the VIMOS (VLT Deep 
Survey) and COMBO 17 classifications. Synthetic models con- 
necting the X-ray and IR SED of AGN indicate that the emis- 
sion arising from pure AGNs should be negligible (i.e., < 10%) 
in high redshift sources detected at 24 /mi ( |Silva et al~| (2004 )). 
To have more precision about the possible fraction of A GN in 
our sample we use a diagnostics colour-colour given by Stern 
et al. | (|2005b and t he diagnostics of the mid-IR (MIR) slope of 



Brand et ai~[ ( |2006| ) (see section 3.3) ; we also use a specific out- 



put of the code CIGALE (see section 4.2). These studies will be 
described later in this paper. 

We decide to distinguish our SED-fitting analysis between 
the sub-sample detected at 70 yum and the whole sample. 
Hereafter we will use ' 70 yum sample' for galaxies detected at 
70 yum and 'total sample' for the whole set of galaxies. 



Web address to use CIGALE : http://www.oamp.fr/cigale/ 



MUSYC collaboration : http://www.astro.yale.edu/MUSYC/ 
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Filters 


Wavelength 


Instrument- Survey 


Galaxies detected 


Mid-IR, IR 


MIPS 1 


24 yi/m 


MlPS/Spitzer - FIDEL 


181 


MIPS 2 


70 jjm 




62 


Near-IR 


IRAC 1 


3.6 jim 


IRAC-GTO Spitzer 


177 


IRAC2 


4.5fim 




161 


IRAC 3 


5.8yL/m 




177 


IRAC 4 


8.0/im 




161 


UV - Optical 


U 


0.35 jim 


MUSYC 


179 


U38 


0.36 iim 




180 


B 


0.46 nm 




181 


V 


0.54 fjm 




181 


R 


0.65 fim 




181 


I 


0.87 fim 




181 


z 


0.90 fim 




181 


J 


1.2/im 




180 


H 


1.6 /urn 




155 


K 


2.1 /urn 




181 


UV 


NUV 


2310 A 


GALEX-DIS 


150 



Table 1. The observed data : wavelengths corresponding to the filters identification of instrument and survey. Our sample contains 
181 galaxies of which 62 are detected at 70 fim. 



SED-fitting techniques 
CIGALE 



description of the code 



The best way to derive physical parameters like the SFH is to fit 
the observed SED with models from a stellar populations synthe- 
sis code. We use the code CIGALE which provides physical in- 
formations about galaxies by fitting the ir UV-to-IR SED. A first 
version of the code was developed by Burgar ella et al. (2005]) 
to reproduce the dust- attenuated stellar emission from the FUV 
to the NIR. Ldust could b e provided as an input parameter only. 
Hence, [Noll et al. (2009b]) especially extended the code to the 
far-IR (FIR) by the consideration of dust emission models to al- 
low a consistent treatment of dust effects on galaxy SEDs. A 
Bayesian-like anal ysis is used to derive galaxies properties and 
Noll et al. (2009b]) applied it to local galaxies. This work repre- 
sents the first attempt to apply CIGALE on galaxies at intermedi- 
ate redshift. CIGALE is based on the use of a UV-optical stellar 
SED plus a dust IR emitting component. First , UV-to-IR models 
are built and secondly, these models are fit to the observed SEDs. 
Practically the code fits the observed data in the UV, optical and 
NIR with models generated with a stellar populations synthesis 
code, assuming a SFH and a dust attenuation scenario as input. 
The energetic balance between dust-enshrouded stellar emission 
and re-emission in the IR is carefully conserved in combining 
the UV/optical and IR SEDs. The IR SEDs are built from the 
Dale & HeloiT|([2002]) templates (hereafter DH02). We refer the 
reader to |Noll et al. | ( |20 09b) for more details. 

In the following we will describe parameters crucial to this 
study. The input parameters values are listed in Table [2]. 



3.1. Stellar populations and star formation histories 

In its current version CIGALE provides two stellar populations 
synthesis models : Maraston et al. (2005) (he reafter M05) and 
PEGASE QFioc & Rocca-Volmerange~| ( [T997] )). Both include a 
treatment of thermally pulsating asymptotic giant branch (TP- 
AGB) stars but in a different way. For instance, in M05 models 
young stellar populations are modelled with the Geneva stel- 



lar evolutionary tracks, while the Padova tracks are used in 
PEGASE. 

TP-AGBs are particular important for a reliable M+ deter- 
mination at high z (M05). The difference between PEGASE and 
M05 models lies in the weight given to the TP-AGBs. Their con- 
tribution will be systematically higher in M05 models. Models 
which do not consider enough TP-AGB need a massive contri- 
bution of old stars in order to match the observed NIR fluxes 
and thus overestimate M+ (M05). In a recent work |Tonini et al. ] 
( 2009]) found that the use of PEGASE does not reproduce IRAC 
luminosities for nearly passive and star-forming galaxies around 
z ~ 2. For these reasons we decide to use the M05 models. 

Two stellar components have been found necessary to recon- 
struct more accurate SFRs in actively star forming galaxies (e.g 
|Erbetal. |p006] ) " 



Lee et al. 



(2009 )). CIGALE uses the single 
stellar populations of M05 and combines two stellar populations 
to reproduce an old and a young stellar population. M05 models 
assume SFH with exponentially decreasing SFR : t\ and T2 rep- 
resent the e-folding time of the old and young stellar population, 
respectively. The SFR is expressed as follow: 



SFRi(0 = SFRo,i.e~ 



SFR 2 (0 = SFR , 2 .e" 



(1) 



(2) 



t\ and t2 represent the age of the old and the young stellar pop- 
ulation and SFRo,i and SFRo,2 are the values of the SFR at t=ti 
and t=t2, respectively. 

We study a sample of LIRGs lying at z=0.7. At this redshift 
the universe was 7 Gyr old with our cosmology, so it makes no 
sense to consider a range up to more than 7 Gyr for the ages of 
both stellar populations. After several tests we realized that we 
cannot make a precise estimate of the e-folding time of the young 
stellar population, so we only consider a SFR constant over the 
time. We decide to fix the age of the old stellar population t\ -1 
Gyr and we suppose 3 possibilities for t\ : T\ = 1,3 or 10 Gyr. 
We choose t<i in the range 0.025 to 2 Gyr and we consider T2 = 20 
Gyr which corresponds to a burst constant over this time-scale. 
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SFR(t) 




old stellar 
population 



population 



Fig. 1. Scenario of SFH : two bursts, one for the old stellar pop- 
ulation having a moderate e-folding time and one for the young 
stellar population constant over time. 



The two SFH components are linked by their mass fraction. 
The fraction of the young stellar population (hereafter f y g p) cor- 
responds to the fraction of the young stellar population mass 
over the total mass and is thus comprised between and 1 . We 
chose values of f y sp in this interval to allow the choice of the 
best fitting value, as we do not have any prescription for this pa- 
rameter. CIGALE gives an instantaneous value of logioSFR as 
an output, defined as log [(1-/J,sp)*SFRi + /J,sp*SFR2]. In the 
following, any occurrence to the SFR will refer to this formula. 

3.2. Dust attenuation and infrared emission 

To modelise the attenuation by dust, the code CIGALE uses at- 
tenuation law of |Calzetti et al. | ([2000 ) as a baseline, and offers 
the possibility to vary the attenuation law and to add a bump. 
We decide to reproduce a pure Calzetti law and to add no bump. 
The stellar population models have to be attenuated before the IR 
emission can be added depending on the dust-absorbed luminos- 
ity. The code allows to consider that young and old stars are not 
affected by dust in the same way. Dust attenuation in the V band, 
Ay, and the dust reduction factor f v are used for the young and 
the old stellar population respectively. The stellar light heats the 
dust which reemits the energy in MIR and FIR. At these wave- 
lengths we have few available data to reconstruct the SED. To 
fit IR observations CIGALE uses semi-empirical one-parameter 
models of DH02 composed of 64 templates. These models are 
parametrised by a, the power law slope of the dust mass over 
heating intensity, defined as follows : 



dM d (U) = U a dU 



(3) 



where Md(U) is the dust mass heated by a radiation field at in- 
tensity U ; a is directly related to f^ m lh^m ratio flux, where 
feo^m and /loo^m represent fluxes of the SED at 60 and 100 yum, 
respectively. 

DH02 provide SEDs corresponding to different exponents a, 
which is a free parameter only for the 70 /im sample. The total 
sample only has 24 jim fluxes, so we take the median value of a 
found by the code for the 70 jum sub-sample. We choose a in the 
interval [1 ; 2.5] to cover a large domain of activity (the activity 
increases as a decreases). 



CD 
< 



CD 



1.0 



0.5 - 



0.0 - 



$ o 



0.5 - O 



■1.0 



o o 




-1 



1 
[5.8-8.0] (AB) 



Stern et al. 



(2005) to 



Fig. 2. Color-color diagram preconized by 
identify AGN candidates. We compare the color [3.6-4.5] versus 
[5.8-8.0] in AB magnitudes for our LIRGs; the area defined by 
the dashed lines represents the region of the color-color diagram 
populated by AGN. 



3.3. Active galactic nuclei contribution 

LIRGS have large Ldust corresponding to an active phase of 
dust enshrouded star formation and/or AGN activity. The ener- 
getic balance between dust-enshrouded stellar emission and re- 
emission in the IR can be violated if dust emission caused by 
an AGN is present. In particular, highly dust-eshrouded AGNs 
represent a problem since they are difficult to identify in the UV 
and optical. Consequently, galaxies with such an AGN contri- 
bution may look like normal star-forming galaxies in the FUV- 
to-NIR wavelength range excepting a warm/hot dust emission in 
the MIR. 

In this paper we are interested in the M* building by LIRGs. 
Therefore it is essential to separate the contribution of AGN and 
starbursts to L^st- The AGN contribution is likely to be domi- 
nant at only very high luminosity (Ldust ^ 2.10 12 L Q , Sanders 
|& Mirabel"] ( | 1996) )). Nevertheless, a small fraction of the total 
IR luminosity can be due to the presence of an AGN even in 
star-formation dominated cases (e.g.|Genzel et al. |fl998 )). 



To disentangle IR dust emission components ca used by stel- 
lar and AGN radiation, we apply the diagnostics of Stern et al. 
(2005 ) to our observations. They combine spectroscopic obser- 
vations with MIR observations of nearly 10000 sources at < z 
< 2 and show that MIR photometry provides a robust technique 
for especially identifying broad-line AGN (90% of broad-line 
AGN and 40% of narrow-line AGN are identified). 

Fig. |2| shows t he colour-colour distribution proposed by 
|Stern et al~| ( |2005| ) in the NIR. The dashed lines represent the 
boundary of the area which characterise AGN features. 15 % of 
the objects are found within this area. However these colou r cri- 
teria may omit AGNs at z « 0.8 an d 2 (|Stern et al. 1(2005) ) . 

We also use the diagnostics of Bran d et al. ( |2006| ) who ex- 
animate the MIR slope £ of the SEDs for their sample. £ corre- 
sponds to the observed flux ratio logio(v/ y (24 yum)/v/ y (8 jim)). 
The MI R slope is supposed t o be steeper for starburst-dominated 
sources. [Brand et al"~| ( |2006| ) assume that sources (z>0.6) with £ 
peaking at « 0.0 and « 0.5 are AGN- and starburst-dominated, 
respectively. Fig.|3]shows the flux ratio £ for the total sample. For 
the vast majority of objects t> is found around 0.5 (-17% satis- 
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Fig. 3. The empty histogram illustrates the flux ratio £, = 
logio(v/ y (24 yum)/y/ y (8 jim)) for the total sample. The black his- 
togram shows the galaxies identified as AGNs by the diagnostic 
of Stern et al. (2005). 
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Fig. 4. SED of the adopted AGN template, with L = 1O 12 L , 
R = 125 pc, and Ay = 32 mag. 



fies £ « 0.0) which lets us think that starbursts dominate the MIR 
emission of our sample. We overplot in black the histogram of 
the objects previously identified as AGNs by Stern et al. (2005). 
We find that they are not necessarily identified by Brand et al. 
(2006). 

CIGALE allows to estimate the fraction of Lg ust due to an 



AGN. Th e code uses AGN templates from |Siebenmorgen et al. 
P004b|a| who provide almost 1500 AGN models differing in 
the luminosity of the non-thermal source, the outer radius of a 
spherical dust cloud covering the AGN, and the amount of atten- 
uation in the visual caused by clouds. All models can be fed into 
CIGALE. Focusing on SEDs providing PAH-free MIR emission, 
the number of suitable models signicantly decreases. As the ref- 
erence model, we take L = 1O 12 L , R = 125 pc, and Ay = 32 
mag (see Fig. [4]). The parameter in our runs which probes the 
contribution of AGN radiation to Ld us t is named /agn for frac- 
tion of AGN . As input we propose values for /agn ranging from 
0.0 to 0.999. The results will be presented in section 4.2. 

Later we will decide to flag objects identified as AGN by at 
least 2 diagnostics ( |Stern et al~| ( |2Q05| ) , |Brand et al~| p006| ), 
CIGALE). 



3.4. Do we retrieve reliable parameters? Construction of a 
mock catalogue dedicated to the sample. 

The code gives an estimate of the basic input parameters plus 
several additionnal output parameters that depend on different 
basic model properties. We mostly focus on SFH parameters for 
the young stellar population. 

We perform a Bayesian analysis to estimat e the physical pa- 
rameters. We use a method similar to the one of Kauffmann et al. 
(2003a). They derive a probability distribution function (PDF) 
for each parameter and estimate the most probable value. We 
use the mean value of this distribution. The expectation value of 
each parameter < x > is given as : 



< x >- 



and the standard deviation cr is derived by : 



(4) 



Ef=i Pi(x~ <x>) 



JtiPi 



(5) 



This method is described as the 'sum' method in INoll et aT~l 
( |2009b| ). 

When we use a numerical code based on SED-fitting, we ob- 
tain results depending on a statistical analysis. The accuracy of 
the estimates is directly dependent on the input fluxes (number 
and quality). We must check which parameter the code is able to 
estimate correctly. Theoretically the Bayesian analysis provides 
an estimation of the reliability of each parameter thanks to the 
shape of their PDF (if the shape is not a Gaussian curve the pa- 
rameter is not well determined ( Wal cher et al. | ( |2008] )), and their 
relative error cr x . 

A more straightforward approach to check the parameter es- 
timation is to generate a mock catalogue from the real sample. 
Our general strategy is to build a specific mock catalogue for 
each set of data (in this paper the LIRGs sample). It is made 
of artificial galaxies for which we do not only know the flux at 
each wavelength but also values of attenuation parameters and 
parameters of the SFH. 

To construct our mock catalogue we follow 3 steps. The first 
step is to run the code on the data to determine the best fitting 
model for each object by a simple^ 2 minimization method. Each 
corresponds to a SED. In the second step we add an error ran- 
domly distributed according to a Gaussian curve with <x=0.1 to 
each flux ; it means that we add an error to the flux which is typ- 
ically 10% of the flux. We have a catalogue of artificial galaxies, 
detected in the same filters as the observed galaxies. Then the 
last step is to run the code on these simulated data and to com- 
pare the exact values of the physical parameters with the values 
estimated by the code. 



4. Results of the SED-fitting analysis 

We perform SED-fitting on our LIRGs sample and on the related 
mock catalogue in order to know which parameter we are able 
to correctly estimate. In the following sections we will describe 
the results of the SED-fitting of the mock catalogue with and 
without the band at 70 yum and then we will discuss results of the 
SED-fitting of the real data. 
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Parameters 


Symbol 


Range 


Attenuation 


V-band attenuation 


Ay 


0.15; 0.30; 0.45; 0.60; 0.75; 0.90; 1.05; 1.20; 1.35; 






1.5; 1.65; 1.8; 1.95; 2.1 


Reduction of A v basic for old SP model 


fv 


0.0; 0.25; 0.50; 0.75; 1.0 


IR SED 


IR power-law slope 


a 


1.0; 1.25; 1.5; 1.75; 2.0; 2.25; 2.5 


AGN-related fraction of L dusU 


Iagn 


0.0; 0.001; 0.005; 0.01; 0.05; 0.1; 0.5; 0.7; 0.999 


Star Formation History 


metallicities (solar metallicity) 


Z 


0.02 


r of old stellar population models in Gyr 


T\ 


0.1; 3.0; 10.0 


ages of old stellar population models in Gyr 


h 


7 


r of young stellar population models in Gyr 


T 2 


20.0 


ages of young stellar population models in Gyr 


t 2 


0.025; 0.05; 0.075; 0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 






0.8; 0.9; 1.0; 1.5; 2.0 


fraction of young stellar population 


fySP 


0.0; 0.001; 0.005; 0.01; 0.05; 0.06; 0.07; 0.08; 0.09; 




0.1; 0.2; 0.3; 0.4; 0.5; 0.6; 0.7; 0.8; 0.999 


IMF 


s 


Salpeter 



Table 2. List of the input parameters of the code CIGALE and their selected range for a specific run. 



4.1. Analysis of the mock catalogue 



We construct two mock catalogues as described in Sect. 3.4 
based on the 70 jim sample (sample detected at 70yum) only. The 
first catalogue is made of galaxies detected at 70 yum ( 70 yum 
mock sample) and the second one is made with the same galaxies 
for which we remove the flux at 70 jim (mock sample) in order 
to know if the fact of having 70 yum flux has an influence on the 
determination of the parameters of the SFH, or M+, L^ust and 
SFR. 

In Fig. [5] is plotted the difference in magnitudes between the 
flux from the best x 2 model estimated by CIGALE and the ob- 
served flux for each broad band filter. Mean and median values 
are nearly equal to zero which means that fluxes are well repro- 
duced by the code to better than 0.1 mag. 

Even if the code is able to well reproduce the observed data, 
we cannot settle for this result. Indeed, it is possible that some 
SEDs are degenerated and correspond to different parameter val- 
ues. We have to compare each input parameter with its Bayesian 
estimation by the code. Fig. [6] and [7] show the input parameters 
of the mock galaxies versus the values estimated by the code. 
Fig. [6] (Fig. [7]) shows results for the 70 yum mock sample (mock 
sample). Values of the Linear Pearson correlation coefficient r 
are gathered in Table [3] 

For the 70 yum mock sample we find very good correlations 
for M star , SFR, L dust , f AGN , and A v with r > 0.9. For t 2 , f y sp and 
fv the correlation is less satisfying, r is equal to 0.34, 0.52, and 
0.64 for the 70 yum mock sample and 0.20, 0.45, and 0.64 for 
the mock sample. The correlation for the parameter a is good 
with r=0.80 for 70 yum mock sample but very bad for the mock 
sample. In the following we will use the mean value of a from 
the 70 yum sample for the objects without detection at 70 yum. For 
/agn r > 0.9 for 70 yum mock sample and 0.76 for the mock sam- 
ple, with relatively small errors. Thus this parameter is correctly 
estimated by the code. 

According to the analysis of the coefficient r and the 
Bayesian method's error, the reliable parameters are : M star , 
Ld U st, SFR, fysp, A v , a (if 70 yum is available) and /agn- 
Consequently in the next section we will only study these pa- 
rameters. 
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Fig. 5. Difference in magnitudes between the flux of the mock 
galaxies and the flux of the best x 2 model for these objects in 
each band. Red crosses and green crosses represent the mean 
and the median value for each wavelength. 



Parameters 


70 fim mock sample 


mock sample 


logio M star 


0.93 


0.92 


log 10 Ldust 


0.93 


0.95 


logio SFR 


0.93 


0.92 


h 


0.34 


0.20 


fySP 


0.52 


0.45 


A v 


0.92 


0.80 


fv 


0.64 


0.64 


a 


0.80 


0.04 


/agn 


0.92 


0.76 



Table 3. Estimation of the linear correlation coefficient of 
Pearson between the exact value and the value estimated by 
CIGALE for each parameter of the mock catalogue. We esti- 
mate this coefficient for the 70 yum mock sample and the mock 
sample. 
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Fig. 6. On the x axis the value for 8 parameters of the mock galaxies in the case of a detection at 70 yum, and on the y axis their 
values estimated by CIGALE. The red line corresponds to the ratio 1:1. 



4.2. Analysis of the LIRGs sample 

We now present the results for the LIRGs sample itself. Fig. [5] 
shows distributions for parameters of interest, calculated with 
the Bayesian analysis described previously for the sub- sample 
detected at 70 yum and for the whole sample. 

The masses found for the 70 yum sample are shifted towards 
larger masses in comparison with the total sample. We observe 
a similiar situation for Ld us t and the SFR; the values are in the 
range [ 10 11 ; 10 12 ] L and [10 ; 92] M .yr _1 respectively with 
mean values larger than ones found for the total sample. This 
shift is expected because of the 70 yum detection limit (2.7 mJy 
at 5<x) : at this wavelength we detect only luminous and mas- 
sive galaxies. The distribution of f y gp is broad with a long tail 
towards high values. Ay lies between 0.5 and 2.1 mag with very 
few objects under 1.0 mag and a quite homogeneous distribution 
between 1.0 and 2.0. 

Galaxies in the total sample have M* between 10 10 and 10 12 
M Q with a peak at 1O 1O - 8 M . We find the SFR between 3 and 92 
M yr _1 with a peak at 23 M yr _1 . For f yS p, and A v we observe 
the same distribution as for the 70 yum sample. For both samples 
/agn is relatively low, between 0.0 and 0.3 with the majority of 
the objects in the interval [0;0.1]. We consider that there is a real 
contamination of Ldust by an AGN when /agn > 15% because a 



contamination lower thanl5% does not significantly modify the 
total IR emission. We must also account for the uncertainty of 
the /agn determination : we only consider as AGN contaminated 
objects with /agn> 15% and fAGN~o~AGN> where (Tagn is the 
standard deviation (see Eq. 5). 

We i dentify 21 objects in the total sample . The diagnostics of 
Stern et al. | ( |20Q5| ) identifies 31 objects, and |Brand etal. | < |20Q6| ) 



9 objects. However, all the objects considered as AGN by the 
co de are not sys tematically identi fied as AGN in the diagnostics 
of |Brand et aT] ( |2006] ) and |Stern et n~| ( [2005] ). We consider ob- 
jects as contaminated by AGN when they are flagged by two of 
the three diagnostics. We identify 15 objects in the total sample 
of which 6 are in the 70 yum sample ; these objects will be flagged 
in the following plots. 

We have also investigated the small subsample of galaxies 
not detected in NUV (31 over 181 galaxies) and did not find any 
statistically significant difference in their main characteristics. 
The largest effect is found for Ay : Ay =1.75 +/- 0.25 mag for 
the non detections in NUV and Ay =1.33+/- 0.34 mag for the 
detections in NUV. A higher dust attenuation for galaxies not 
detected in NUV is indeed expected. 
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Fig. 7. Same plot as Fig.^in the case of non-detection at 70 fim. 



5. Discussion 

5.1. Luminosity-colour distribution and the implication for the 
dust temperature 

SED shape of a LIRG in the 8-1000 yum range is due to emis- 
sion from dust. The SED rest-frame peak defines the wavelength 
of maximum energetic output (in vf v ) ; it is usually in the 40- 
140 /mi range and the SED is often approximated by a modified 
black-body (B v v^) of T ~ 20-60 K and emissivity 0-1-2. 

Local IRAS (Infrared Astronomical Satellite) galaxies ex- 
hibit correlations b etween the color f^ m lf\^ m and L dus t (e.g 
Dale et al7] ( |2001| )). fm^mlfm^m ratio is often taken as an indi- 
cator of the typical heating conditions in the interstellar medium 
and is indicative of a characteristic dust temperature (hereafter 
T dust) , feofim/fioojim increasing with T dust . 

In the nearby universe LIRGs are characterized by emission 
from dust at high temperatures (Sanders & Mirabel (1996)). 
Fitting single T dust to feoum/fmoum ratio leads to T dust ~20K fo r 
low redshift spirals ( [Reach et al. ] ( [T995] ); [Dunne & Eales |pQQT] » 
and 30 -60K for the h igh-lumin osity o bj ects typically detected by 
IRAS flSoifer & Neugebauerj ; [Stanford et"al~| pOOO] )). 

High-redshift, hyperluminous galaxies can show T dust of up to 
1 10K ( [Lewis et al. HI998) ). 



To reproduce number counts at FIR wavelengths, |Lagache et 
al~[ j2005) , [Chapin et al. | ( |2009a| ), and |Valiante et al. | PQQ9J 
found that it is necessary to take into account cold luminous 
galaxies in the nearby universe. Cold galaxies are found to domi- 
nate the z=0 luminosity function at low luminosities and become 
less important at higher redshifts, because of their passive evolu- 
tion. The fraction of cold galaxies can contribute to the number 
counts up to 50% at 170 /mi. At moderate and high redshift the 
possible impor tance of cold , luminous g alaxi es has also been 
point ed out by Eales et al. ( 1999,2000) and Chapm an et al. 
(2002]). 

We use the fluxes ratio fi^mlfi^m (corresponding to the 
highest wavelengths in our sample) in order to estimate the 
dust temperature. The observed fi\ m lfi^m corresponds to 
funmlUium rest-frame. All our discussions are made in the 
framework of the DH02 models in the [1.0-2.5] a regime (see. 
Table[2]). fu/um/f^m evolves in the opposite way of fe^mlfi^m 
(i.e. Tdust, see Fig. [9]). 

Fig. [TO] exhibits the observed fiA^mlfi^m ratio as a function 

of 

Ldust obtained for the LIRGs in the 70 /im sample. Predictions 
of DH02 models are overplotted. They are calibrated with the 
Marcillac et al. (2006) local relation between f^ m lf\w m (i.e. 
a) and L dus t • For a given L dus t we find that our galaxies have 
a fujimlfAi^im ratio higher than predicted by the models. This 
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Fig. 8. Bayesian results of the code for the following parameters : M star , L^ust, SFR, A v , f yS p and /agn- The empty histogram 
represents the total sample and the full one represents the 70 jim sample. 



implies that they have a lower f^mlfi^m ratio and thus are 
colder than expected. We cannot attribute this result to the pres- 
ence of AGN contaminated objects regarding the small fraction 
of flagged sources. 

We can also directly use the feo^mi 'fioo^m ratio deduced from 
the a parameter estimated by CIGALE. In Fig. [TT] we com- 
pare Ldust and the derived /60/W TioOyum values of our 70 fim 
sample. It will have no sense to make this comparison for the 
rest of the sample for which we estimated a with the median 
value of a found for the sub-samp le. We add to this plot the lo- 
cal Ldust-fauJf i oom re lations of |Chapman et al. | ( [2003] ) and 
Marcillac et al. | ( |2006|), and the lum inosity-temperature mod 
els pro posed by CE01, Chapin et al. | ( |2009a| ), and Valian te et al. 
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shows that 
Z/ioo/zm ver- 



( |2QQ9| ) in the luminosity range [10 11 ; KFTL . Fig 
the code gives a different trend for the relation f^ Ml 
sus Ld U st than both local relations and models. Galaxies having 
the lowest IR luminosity are in agreement with them, but more 
luminous objects are found colder than predicted, and globally 
they evolve in t he opp osite way to local relations and models. 
Cha nial et al. | ( |2007| ) make a selection at 60 and 170 yum us- 
ing respectively IRAS and ISO observations at z=0 and find 
that galaxies selected at 170 yum are colder than ones detected 
at shorter wavelengths. They conclude that selecting galaxies 



at long wavelengths introduce a bias towards cold galaxies. 
Selecting the sample at a wavelength long wards of the peak, by 
definition, does no t predispose towards warm sources ( Sanders 
& Mirabel (1996 )). Our sample is selected by rest- frame -41 
fim emission, i.e. on warm dust; yet, we mostly find colder T^ust 
at a given luminosity than predicted by models of DH02 and 
CE01, calibrated in the nearby universe. This tends to support 
the interpretation that the offset found towards colder tempera- 
tures at a g iven luminosity is at least in part a real difference. 
Zheng et al. (2007 ) suggest that this tendency towards colder 



Tdust reflects a difference in dust and star formation geometry: 
whereas local LIRGs tend to be interacting systems with rela- 
tively compa ct very intense star formation and of comparable 
masses (e.g., |Wang et al. | p006| )), 0.5<z<l LIR Gs tend to be 
disk-d o minated, relatively undisturbed galaxies ([Zheng et al. 
p004] ); |Bell et al. | ( [2005] ); [Melbourne eTal~| ( |2005| )). They sug- 
gest that these disk-galaxies host widespread intense star forma- 
tion, like star formation in local spirals but scaled up, leading to 
relatively cold Tdust- 

Our results are in agreement with[Symeonidi s et al. ( 2009 ). 
They compare FIR properties of LIRGs and Ultra Luminous 
IfraRed Galaxies (ULIRGs, L dust > 1O 12 L ) in the 0.1<z<2 red- 



shift range, based on MIPS and IRAC data with L dus t > 10 iU L ( 
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Fig. 10. This figure illustrates the ratio between the observed fluxes at 24 jam and 70 jim as a function of the L dust estimated by the 



code. The solid line corresponds to the results given by the DH02 templates at z=0.7. The dashed cyan line is given by Siebenmorgen 



et al. (2004b) for typical AGN contaminated object. AGNs described in section 4.2 are shown as red stars. 
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Fig. 9. Variations of fu m IU\ m versus fa m lfvw m (expressed 
in Jansky) for the DH02 templates used in the paper, correspond- 
ing to the a range [1.0 - 2.5]. 



70 yum selected objects, with those of local sources of equivalent 
luminosity from the local IRAS Bright Galaxy Sample. They fit 
their data with the Siebenmorgen & Kriigel ( 2007) theoretical 



SED templates and find that their distant LIRGs and ULIRGs 
are colder than their local equivalents. They found a significant 
cold-dust associated FIR excess not detected in local sources of 
comparable luminosity directly linked to the evolution in dust 
and star formation properties from the local to the high redshift 
Universe. 

Herschel is now providing longward 70 jim data allowing a 
reliable estimation of Tdust- 



5.2. Star formation rate and stellar mass 

How and when galaxies build up their M+ is still a major ques- 
tion in observational cosmology. While a general consensus has 
been reached in rec ent years on the evolution of the galaxy stella r 
mass function (e.g., |Dickinson et aL"| (|2003'); Dr ory et aI7|(|2004| ); 



Bundy et al.| ( |2005] ); ]Pannella |p006| ); |Marchesini etaL| (2008 ) 
the redshift evolution of the SFR as a function of M+ still re- 
mains unclear. The evolution of the galaxy SFR-M* relationship 
provides key constraints on M+ assembly histories of galaxies. 
For star-forming galaxies a correlation is observed between SFR 
and M+ with a slope near to unity from z ~ to 2 and a scatter 
of 3 dex dDaddi et al. | { |2007) ; |Elbaz eTaT| ( [2007l ); |Noeske et| 



al 



(2007a)). The slope and scatter are expected to be a direct 
reflection of the similarity in the shape of galaxy SFHs among 
various m odels and at various masses ( [Dave | ( |2008| ), |Noeske et| 
|al. | ( |2007b| )). It is worth noting that other z = 2 observations sug- 
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Fig. 11. Comparaison of the relation log R(60/100) - logioL dust given by IR SEDs of DH02 after the fit to the data, with the values 



et al. 


(2006 


et al. 


(2009 



, respectively. We overplot models of CE01 (green dotted line), Chapi n~et al. | ( |2009a| ) (dashed black line) and Valiante 
(dash-dotted black line). AGN detections are shown as red stars. 



gest a somewhat large scatter ( Shap ley et al. | ( [2005] );[Papovich 
et al. (2006 )), so the observational picture is not entirely settled. 



Fig[T2|shows the SFR of each galaxy of our sample as a func- 
tion of the corresponding M*, and the results of various models 
and observations : analytical models of Buat et aL] (IT008 ) and 



Noeske et al. ( 2007a|), the sem i-analytical Millennium simula- 
tions of |Kitzbichler & White |p007) used in|Buat et al.|p008) 
and observations of|Elbaz et al. | ( |2007] ), |Daddi et al. | ( |2007| ),"and 



Santini et al. 



(2009 



Our sample of LIRGs has an average red- 
shift of ~ 0.7. We find the vast majority of the sample between 
the correlations at z=0 and z=l given by Elba z et aL~| ( |2QQ7j ). Our 
data cover the extreme part of the Millennium simulations where 
we can find galaxies with the largest SFRs between ~ 10 - 100 
M yr _1 . We do not observe a correlation between SFR and M+, 
but instead a flat distribution of SFR. We consider a tight range 
of SFR since we select objects with L^ust >10 n L© (see below) 
: this range is not sufficient to search for a general correlation 
between M* and the SFR. However in the considered SFR range 
we observe a large scatter in M+ . Indeed galaxies whose young 
stellar population is younger than 300 Myr exhibit a larger dis- 
persion in M+ than galaxies whose SFR is constant over a longer 



period (Fig 12). Therefore the SFH of our objects has an influ- 
ence on the dispersion of the SFR-M* plot. This is in agreement 
with the conclusions of |Dave | ( |2008| ) and |Noeske et al. | ( |2007b| ). 



The specific star formation rate (SSFR) defined as the ratio 
between SFR and M+ (e.g., Kennicutt et al. (2005)) is also com- 
monly used to trace the SFH in galaxies. The SSFR is plotted in 
Fig[T5]as a function of the corresponding M+ and we can see the 
increase of the SSFRs towards lower M*. Before any interpre- 
tation we can estimate a detection limit in SSFR : L^ust — 10 u 
L Q translates into SFR //m = 17.78 M yr _1 at z = 0.7 ( Eq. (1) 
of Buat et al. (2008)). We overplay this detection limit and see 
that the trend of SSFR with M* is largely constrained by this 
limit. We compare our results with the same models and obser- 
vations as illustrated in Fig 12 SSFRs of our intermediate-mass 
galaxies (M ir <lO n M Q ) are in agreement with those observed by 
Sant ini et aI7| ( |2009| ), but for these values of M* we cannot detect 
more quiescent galaxies because of the detection limit. SSFRs of 
large-mas s galaxies (M*> 10 n M o ) are better described by the 
models of |Buat e t al. (2008). These massive galaxies have SSFRs 



which do not exceed 2.5 10" iU yr _1 (corresponding to a SFR-25 
M©yr _1 ) and never reach the SSFRs of intermediate-mass galax- 
ies observed in this study. In order to distinguish galaxies for 
which the young stellar population significantly contributes to 
M*, we identify objects for which f y $p estimated with CIGALE 
is larger or lower than 10% and see that the 2 groups are well- 
separated according to their stellar mass: galaxies with M+ < 
10 n M o have more than 10% of their mass due to the young stel- 
lar population, whereas in galaxies with M* > 10 u M o this mass 
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Fig. 12. The relation SFR-M* for the total sample. Filled circles (empty circles) illustrate galaxies for which the age of the young 
stellar population is found greater (lower or equal) than 0.3 Gyr by CIGALE. AG Ns are represented by red stars. We overplot the 
semi- analytical Millennium simulations (cyan diamonds), analytical models of Buat et al. ( 2008) (dash-dotted yellow line) at z=0.7 



and Noeske et al. 



(12007 a| ) (connected blue crosses) at z~0.7, and observations of Elba z et al. | ( |2007] ) (solid black line, z=0 and 1), 
paddi et al.](|2007|l (solid black line z=2) and|Santini et al. |(|2009l) at z ~ - 7 (dashed green line). 



represents less than 10%. This result is in agreement with the 
common picture that massive galaxies formed the bulk of their 
stars early and on shorter timescales, while numerous less mas- 
sive galaxies evolve on longer time-scales, a phenomenon gener- 
ally linked to the downsizing s cenario (e.g. Cowie et al. (1996), 
Brin chmann & Ellis |(|2000l) ;|Fontana et al | ( [2003]) ; |Feulner 
et al. | ( |2005|) ; |Perez-Gonzalez et al. | ( |2005| ) ; |Papovich et"aL 



p006| ) - pamen et al. | ( |2009] )): 



6. Conclusions 

We performed a SED-fitting analysis of a sample of LIRGs at 
z=0.7 detected at 24 jum by S pitzer/MlPS and for which 80 % of 
the objects have a detection at 231 nm by GALEX. This sample 
is observed over a large range of wavelengths, from NUV to FIR 
(70 yum) and is made of 181 galaxies which 62 are detected at 
70 fim. We fit the SEDs of our galaxies with the CIGALE code 
( [Burgarella eTaE] ( [2005] ), |Noll et al. | ( |2009b| )) which combines 
stellar and dust emissions in a physical way. This study is the first 
use of CIGALE at redshift larger than 0. The stellar populations 
synthesis code of Maras ton et al. | ( |2005| ) is adopted to model 
the stellar emission (UV, optical and NIR wavelength s). Then 
the code uses semi-empirical one-parameter models of D ale &| 
Helou (2002) to reproduce the dust emission in MIR-to-FIR 



CIGALE allows the estimation of several physical parame- 
ters based on a Bayesian-like analysis. We show that we are able 
to estimate reliable stellar ma sses, star formation r ates, fractions 
of young stellar population, Dale & Helou ( 2002) IR-templates 
parameter a (if 70 fim is available), dust attenuation, dust lumi- 
nosities as well as the fraction of dust emission due to an ac- 
tive galactic nucleus. The investigation is performed by building 
mock catalogues based on the real samples. 

The 70 fim sample exhibits colder dust temperature s (as 
traced by the rati o L24 fim /Ljo fim and the a parameter of the Dale 
& Helou ] ( |2002| ) templates) than expected from local relations 



wavelengths. 



between dust luminosity and temperature, confirming other re- 
cent results for similar samples. 

Our LIRGs appear to form stars actively; they exhibit a flat 
distribution with a large scatter in the star formation rate - stellar 
mass plot. The amplitude of the dispersion is related to the age 
of the young stellar population, a tighter distribution being found 
for the largest ages. 

We find that our galaxies with a stellar mass > lO n M have 
less than 10% of their mass coming from the young stellar pop- 
ulation. The specific star formation rate for these massive galax- 
ies never reaches the one found for intermediate-mass galaxies 
of our sample, confirming the downsizing scenario. 

The multiwavelengths data analysis performed in this study 
provide reliable estimates of several physical parameters but 
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Fig. 13. The SSFR-M* relation for the total sample. Green dots and yellow dots represent objects for which f y s p are larger or lower 
than 10%, respectively. AGNs contaminated objects are represented by red stars. Black dashed line shows the detection limit. We 

]J (|2007a| ) (connected blue 

crosses) at z~ 



overplot analytical models of |Buat et al.|(|2008|) (dash-do tted yellow line) at z=0.7 and Noes ke et al. 

( |2007| ) (solid black line, z=l) and |Santini et al. | ( |2009t at z ~°- 7 (dashed green 



line). 



-0.7, and observations of |Elbaz et al. 



may turn out insufficient to determine accurate dust tempera- 
tures. Forthcoming data from Herschel will help us to better 
constrain galaxies SEDs and thus to derive more reliable param- 
eters. 
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